Observation of dipole-mode vector solitons 
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We report on the first experimental observation of a novel type of optical vector soliton, a 
dipole-mode soliton, recently predicted theoretically. We show that these vector solitons can be 
generated in a photorefractive medium employing two different processes: a phase imprinting, and 
a symmetry-breaking instability of a vortex-mode vector soliton. The experimental results display 
remarkable agreement with the theory, and confirm the robust nature of these radially asymmetric 
two-component solitary waves. 



Optical spatial solitons in (2+1) dimensions are 
particle-like solitary waves propagating in a nonlinear 
bulk medium pi. The exhaustive research of the past 
decade has shown that these "light particles" can pos- 
sess topological phase properties analogous to a charge. 
Moreover, several light beams can combine to produce 
a vector soliton with a complex internal structure. This 
process can be thought of as the formation of a "solitonic 
molecule" from the constitutuents of different charge. 

Recently, the existence of the most robust "solitonic 
molecule" , the dipole-mode vector soliton, has been pre- 
dicted B. This novel optical vector soliton originates 
from trapping of a dipole-mode beam by a waveguide 
created by a fundamental soliton in the co-propagating, 
incoherently coupled, beam. While many other topolog- 
ically complex structures may be created, it is only the 
dipole mode that is expected to generate a family of dy- 
namically robust vector solitons [0. The closest counter- 
example is the vortex-mode vector soliton Q which has a 
node-less shape in one component and a ring-like vortex 
in the other component. This radially symmetric, vector 
soliton undergoes a nontrivial symmetry-breaking insta- 
bility which transforms it into a more stable object 
- radially asymmetric dipole-mode vector soliton, even in 
an isotropic nonlinear medium. 

While the existence and robustness of the dipole-mode 
vector solitons have been established theoretically for a 
general model of an isotropic medium with saturable non- 
linearity the main question still stands: Is the stabil- 
ity of these asymmetric solitons, as opposed to the radi- 
ally symmetric vortex-mode solitons, a fundamental phe- 
nomenon that can be demonstrated experimentally? 

In this Letter we answer this question positively. We 
observe dipole-mode solitons in strontium barium nio- 
bate (SBN) photorefractive crystals experimentally, by 
employing two different techniques. First, we use a spe- 
cially fabricated phase mask to create a dipole- like struc- 
ture in one of the co- propagating, mutually incoherent, 
beams. Second, we observe the symmetry breaking of the 
vortex-mode soliton and the formation of a dipole-mode 
soliton, as predicted by the theory. 

Theoretical results. We consider two incoher- 



ently interacting optical beams propagating in a bulk, 
isotropic, saturable medium. The model describes (2+1)- 
dimensional screening solitons in photorefractive (PR) 
materials in the isotropic approximation j^]. It repre- 
sents a great simplification with respect to a more realis- 
tic treatment taking into account the inherent anisotropy 
of the nonlocal nonlinear response of a PR crystal . 
However, we will show that it does provide correct quali- 
tative predictions of the phenomena observed experimen- 
tally. Moreover, due to the generality of the model, one 
can expect to observe similar phenomena in other non- 
linear isotropic or weakly anisotropic systems. 

The normalized equations for the slowly varying beam 
envelopes, E\ and E2, can be written as follows ^,|[: 
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where Aj_ is the transverse Laplacian, and z is the prop- 
agation direction. Stationary, (2+l)-dimensional solu- 
tions of these equations can be found in the form E\ = 
47h u(x, y) exp(i^iz), E 2 = y/Jh w(x, y) e-xp(i[3 2 z), 
where (3\ and [3 2 are two independent propagation con- 
stants. Measuring the coordinates x and y in the units 
of VPi) an d introducing the soliton parameter A — 
(1 — /?2)/(l — Pi), we derive the stationary equations for 
the normalized envelopes u and w: 



A ± u-u + F(I)u = 0, 
A ± w - \w + F(I)w = 0, 



(2) 



I = u 2 + w 2 , 
The limit s 



and s = 1 — Pi is 
— > corresponds 



where F(I) = 1(1 + si)- 1 
the saturation parameter 
to the Kerr medium [^J. 

As has been recently shown in 0, when one of the 
beams components, say w, is weak, a soliton formed by 
the it-component induces a change in the refractive index 
of the PR material that traps and guides the weaker w- 
component. If the waveguide is induced by a fundamen- 
tal, bell-shaped soliton, its nonlinear guided modes form 
a set analogous to the Hermite- Gaussian and Laguerre- 
Gaussian linear modes supported by a radially symmet- 
ric waveguide ||. At higher intensities of the trapped 
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beam, the two beams form a vector soliton that is self- 
trapped by a composite refractive index change induced 
by both beams. Apart from the bell-shaped vector soli- 
tons generated by the ground-state mode of the soliton- 
induced waveguide, one can anticipate various types of 
higher-order (2+1) -dimensional vector solitons. For ex- 
ample, a radially-symmetric vortex-mode vector soliton 
HQ is generated by a node-less beam guiding the compo- 
nent with a nonzero topological charge and the doughnut 
structure of a Laguerre-Gaussian (LGj) mode. A radi- 
ally asymmetric dipole-mode vector soliton has one 
node-less component, while the other component has a 
structure of a Hermite-Gaussian (HGqi) mode. 









FIG. 1. Examples of constituents of (a) dipole-mode and 
(b) vortex-mode vector solitons for s = 0.3 and s = 0.65, re- 
spectively (top and middle row) shown with the phase distri- 
butions needed to generate the w-component experimentally 
(bottom row). 



The existence domain for the radially asymmetric 
dipole-mode vector solitons of our model (0) was deter- 
mined in Ref. 0. A typical example of such a soliton is 
shown in Fig. 3(a), along with a phase distribution in a 
dipole-mode constituent. For any given s, the solutions 
are characterized by a certain cut-off value A, when the 
dipole-mode component vanishes. Near the cut-off, the 
vector soliton can be approximately described by the lin- 
ear waveguiding theory With increasing A, the w compo- 
nent grows and deforms the effective waveguide generated 
by the u-mode, so that the u-componcnt elongates and 
becomes radially asymmetric. The linear and dynamical 
stability analysis of these solitons has revealed their as- 
tounding dynamical stability, with respect to both small- 
and large-amplitude perturbations p| . Overall, the anal- 
ysis performed in suggested that the dipole-mode vor- 
tex solitons would lend themselves to experimental ob- 
servation more easily than their unstable vortex-mode 



counterparts. 

The family of radially symmetric vortex-mode solutions 
of Eqs. (^J), in the form: u = u{r), w = w(r) exp(imcj)) , 
where r — y/ (x 2 + y 2 ), has been numerically found in 
@,||- The u component of this solution has no topologi- 
cal charge and the w component carries a single-charged 
vortex (m = ±1). In Fig. |l|(b), we present a typical so- 
lution of this family, along with the helical phase distri- 
bution needed to generate a vortex in the iu-component. 
The dynamical and linear stability analysis of these solu- 
tions conducted in ^,|| has proved that all such vortex- 
mode vector solitons are linearly unstable. The instability 
growth rate is positive for any vanishingly small ampli- 
tude of the w component, and it increases rapidly with 
the vortex intensity. 

Remarkably, an unstable vortex-mode soliton displays 
a symmetry-breaking instability. It always decays into a 
radially asymmetric dipole-mode soliton with nonzero an- 
gular momentum which can survive for very large prop- 
agation distances 0. The break-up of the vortex-mode 
vector soliton in a saturable medium is a nontrivial ef- 
fect which, as we show below, can be used to observe 
dipole-mode vector solitons in PR crystals. 
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FIG. 2. Experimental setup. 

Experimental techniques. Formation of dipole mode 
vector solitons has been investigated using an experimen- 
tal setup shown schematically in Fig. |^. An extraordi- 
nary polarized laser beam (A = 532nm) was split into 
two parts. First beam was then transmitted through 
the phase mask (or glass slide) in order to imprint the 
required phase structure. In this way, we could obtain 
either an optical vortex with intensity vanishing in the 
center of the beam [see Fig. 0(b)] or a dipole-like struc- 
ture with a phase jump of n across the beam along its 
transverse direction [see Fig. 0(a)] that is perpendicular 
to the optical axis of the crystal [the axis (c) in Fig. 0. 
The second beam was transmitted through the system of 
spherical or/and cylindrical lenses in order to form cither 
a circular or elliptically shaped spot of desired size. The 
beams were later combined using the beam splitter BS2 
and focused into the input face of the photorefractive 
crystal. We used two samples of the ferroelectric SBN 
crystal doped with Cerium (0.002% by weight). Their 
dimensions (ax bx c) are (10 x 6 x 5) mm or (15 x 8 x 5) 
mm. It is well known that photorefractive crystals bi- 
ased with strong DC electric field exhibit strong positive 
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or negative nonlinearity depending on the polarity of the 
field pof . In our case, the SBN crystal was biased with 
the DC field of 1.5 — 2.5 kV applied along an optical axis 
[the axis (c) in Fig. 0]. The resulting photorefractive 
nonlinearity was of a self-focusing, saturable character. 
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FIG. 3. (a) Results of numerical simulations showing 
the generation of a stable dipole-mode soliton (s = 0.65, 
A = 0.5, z = 35). (b) Experimental results showing formation 
of a dipole-mode soliton (bottom) from the two out-of-phase 
beams (top, right) and co-propagating Gaussian beam (top, 
left). Experimental parameters are: V = 2 kV , z — 15 
mm, and the initial powers of the u and w components 
Pu = P w = 0.6 nW. 



made both beams effectively incoherent inside the pho- 
torefractive crystal as the slow photorefractive response 
does not follow fast changes of the relative phase between 
the beams. In all our experiments the initial power of the 
beams did not exceed a few microwatts. 
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FIG. 5. Suppression of the vortex filamentation due to 
modulational instability in the presence of a co-propagating 
fundamental- mode beam (see discussion in the text). Numer- 
ical results are obtained for s = 0.65 and A = 0.5: (a) z = 35, 
(b) z = 90. 
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FIG. 4. (a) Results of numerical simulations showing 
strong repulsion of the input dipole-mode lobes without the 
second component, (b) Experimental results showing the 
strong repulsion the two out-of-phase solitons formed by 
the dipole lobes without the co-propagating Gaussian beam. 
[Theoretical and experimental parameters as in Fig. ^(a,b)]. 

To control the degree of saturation, we illuminated 
the crystal with a wide beam derived from a white light 
source. We estimated the initial (i.e. at the input face 
of the crystal) degree of saturation to be of the order of 
unity in all our experiments. Since both the components 
forming a vector soliton have to be mutually incoher- 
ent, one of the constituent beams was reflected from a 
vibrating mirror mounted on a piezoelectric transducer 
(M2, PZT). Driving this transducer with an AC signal 
of 1 kHz imposed a frequency shift onto the beam. This 




FIG. 6. Experimental demonstration of suppression of the 
vortex filamentation in the presence of co-propagating fun- 
damental-mode beam. Shown are (a) the vortex-bearing 
input beam, (b) output beam after propagating without 
the Gaussian beam, (c) output intensity of the dipole-mode 
soliton formed by two co-propagating vortex and Gaussian 
beams, and (d) intensity of the dipole-mode constituent of 
the two-component vector soliton. Parameters are: V = 2.3 
kV, P u — Pw = 0.3 )J,W, z = 6 mm. 



Experiment vs theory. First, we generate the dipole- 
mode soliton by applying a phase mask to one of the 
input beams to create a dipole-like structure as in Fig. 
0(a). If the dipole-mode bearing beam is launched si- 
multaneously with a Gaussian beam shown in Fig. ||(b) 
(top row), and the two beams are made mutually inco- 
herent, a robust dipole-mode vector soliton is generated 
[Fig. |(b), bottom row]. However, when a dipole- mode- 
bearing beam is launched in the absence of a Gaussian 
beam, two out-of-phase lobes of the dipole beam form two 
coherently interacting fundamental solitons that strongly 
repel each other [Fig. |](b)]. 

Our numerical simulations of the original equations (Q) 
provide excellent agreement with the experiment, as is 
seen in Figs. ||(a) and^(a). The propagation distance of 
z = 35 used in these simulations corresponds to w 9 mm 
propagation distance in the PR crystal, which is a pri- 
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on smaller than the crystal length in most experimental 
runs. 

Next, we examine experimentally a non-trivial break- 
up of the vortex mode in the presence of the Gaussian 
beam. Without the Gaussian beam, a scalar vortex ex- 
hibits filamentation due to the modulational instability 
plE2l. An intermediate stage of a complex vortex fil- 
amentation is shown in Fig. |5|(a). The corresponding 
results of an experimental propagation of a single-charge 
vortex is shown in Figs. ^(a,b). The filamentation pic- 
ture is different in the experimental situation owing to the 
and the inherent anisotropy of the nonlocal nonlinear re- 
sponse of a PR crystal. When the vortex component w is 
launched simultaneously with the uncharged component 
it, the scenario of the instability development changes 
dramatically. The vortex-mode instability leads to the 
formation of a dipole-mode vector soliton, as we demon- 
strate numerically [see Fig. ||(b)] and experimentally [see 
Figs. ^|(c,d)]. Importantly, the phase distribution across 
the profile of the resulting dipole-mode beam calculated 
numerically and measured in the interference experiment 
[see Figs. [?] and 0), confirms that the observed localized 
two-component structure is indeed the dipole-mode vec- 
tor soliton predicted theoretically. 




FIG. 7. Numerical results showing the intensity and phase 
distribution for the dipole-mode component of a vector soli- 
ton formed via decay of a vortex-mode vector soliton (s — 0.3, 
A = 0.6, z = 30). 
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FIG. 8. Experimental results showing (a,c) the intensities 
and (b,d) interference pattern for the dipole-mode lobes gen- 
erated (a,b) via the phase-imprinting experiment, as in Fig. 
^|(b), and (c,d) via the symmetry-breaking of a vortex-mode 
soliton, as in Fig. [j|(c,d). The shift seen in the interference 
fringes of the dipole-mode lobes overlapped with a mutually 
coherent plane wave, (b,d), is the result of the n phase differ- 
ence between the lobes. Parameters are: (a-b) V = 1.8 kV, 
Pu = Pw = 1 1-iW, z = 15 mm; (c-d) corresponds to V = 1.8 
kV, P u = 1.7 fxW, P w = 0.4 fiW, z = 10 mm. 



In conclusion, we have generated experimentally a 
novel type of vector optical soliton in a bulk medium. 
This soliton has a radially asymmetric structure and orig- 
inates from trapping of a dipole mode by the soliton- 
induced waveguide, being much more robust than the 
corresponding vortex-mode vector soliton. Our experi- 
mental results are verified by a systematic comparison 
with the theory. 
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